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I. INTRODUCTION 



The description of hadron form factors within QCD represents a major challenge because 
in such exclusive processes intact hadrons appear in the initial and final states, bearing the 
quark-gluon binding effects controlled by nonperturbative dynamics. Essential to imple- 
menting perturbation theory is the factorization assumption of a short-distance part from a 
large-distance remainder, the latter being mainly nonperturbative. It was proven long ago 
[1, 2, 3, 4] that the pion's electromagnetic form factor can be factorized at large Q 2 in the 
form of a convolution of a hard-scattering amplitude Th and two distribution amplitudes 
(DA) (fx [5] to describe the initial and final pion states. 

Following this type of approach, it has been shown by many authors (see, e.g., [6, 7, 8, 9, 
10] and references cited therein) that the factorizable (hard) part of the electromagnetic pion 
form factor is too small to afford for a good agreement with the available experimental data 
[11, 12]. Admittedly, the quality of the existing data at higher Q 2 values is too poor to draw 
any definitive conclusions, though — at least the tendency — is clear: if at al, the perturbative 
regime seems to be far outside the currently probed momentum scales. For that reason, it 
seems reasonable to pursue alternative methods to compute the pion's electromagnetic form 
factor which do not rely on perturbation theory. 

Our present work is in the context of a three-point AAV (A for axial, V for vector cur- 
rent) correlator, having recourse to the operator product expansion (OPE) and a dispersion 
relation. Our techniques are closely related to the QCD sum-rule approach of [13, 14], but 
with the important difference that we use exclusively nonlocal condensates. In particular, 
a quark-gluon-antiquark nonlocal condensate is employed in which all three inter-parton 
separations are nonlocal. 

Moreover, we use a spectral density which includes terms of 0(a s ). The influence of 
this next-to-leading-order (NLO) contribution to the pion form factor turns out to be quite 
important, reaching the level of 20%. To describe the momentum distribution of vacuum 
quarks, we go beyond the minimal Gaussian model, used before [15, 16, 17], and consider also 
a model, obtained more recently by two of us [18], that has the following methodological 
advantages: (a) it satisfies the QCD equations of motion and (b) it minimizes the non- 
transversality of the VV correlator. 

The principal results of this paper are the following: 

(i) We give a general formalism for the calculation of the pion's electromagnetic form fac- 
tor that contains several methodological improvements relative to all previous approaches 
employing QCD sum rules with vacuum (local and nonlocal) condensates. 

(ii) We suggest a way how to determine the threshold in the local-duality approach for inter- 
mediate values of Q 2 where measurements have already been carried out or are planned. 

(iii) We present predictions for the pion form factor, including also inherent theoretical uncer- 
tainties, that are in good agreement with new lattice results and real experimental data from 
the Cornell and JLab Collaborations. 

The plan of this paper is as follows. In the next section we recall the standard QCD 
description of the pion's electromagnetic form factor according to the factorization theo- 
rem. A brief report on the current status of the obtained predictions is given, avoiding 
technicalities. In Sec. Ill, we discuss the method of QCD sum rules with local and nonlocal 
condensates in the calculation of the pion form factor. The analysis of the present work 
is detailed in Sec. IV. We present our computations and the obtained results and compare 
them with other theoretical predictions, recent lattice simulations, and experimental data. 
This section contains also estimates of the inherent theoretical uncertainties of the applied 
method. In Sec. V we state our conclusions and summarize our main results. 
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II. CONVOLUTION SCHEME FOR THE PION FORM FACTOR IN QCD 



Applying the factorization theorem [2, 4], the pion's electromagnetic form factor can be 
written in the form (for reviews see, for instance, [19, 20, 21]) 



i i 

2 



F AQ ) = // ^(zjA* ) T n(x,y,Q ,fi )(p*(y,[i )dxdy. (2.1) 







The hard amplitude is the sum of all Feynman diagrams in which the struck quark is con- 
nected to the spectator via highly off-shell gluon propagators, meaning that the transverse 
interquark distance is rather small, i.e., of the order of the inverse large momentum transfer 
Q, and that both partons share comparable fractions of longitudinal momentum Xj = pf / P + , 
where the parton four-momentum in light-cone coordinates is = (pf ,p^ , k_u) and P is 
the initial four-momentum of the pion. Momentum conservation under the assumption of 
exact u — d symmetry implies that Y^=i x i = 1 an d Y^f=i kjj = 0. In our case, x\ = x and 
X2 = 1 — x = x. 

The unknown binding effects of the pion state have been absorbed into the pion DA which 
at the leading-twist level two is represented by the valence-state wave function on the light 
cone averaged over transverse momenta up to the factorization scale \i of the process. It is 
defined by the following universal operator matrix element (see, e.g., [19] for a review) 



(0 | d(z)^~f 5 C(z,0)u(0) | tt(P)) 
with the normalization 



= i 

z 2 =0 



f v P" [ dxe**W ^{x,fi 2 ) (2.2) 
Jo 



i 

2^ 



<p v (x,fi Q ) dx = 1 , (2.3) 
where = 130.7 ± 0.4 MeV [22] is the pion decay constant defined by 

(0|J(0)7 M 7 5 n(0)|7r+(P)> = iP,U . (2.4) 

Here 



C(z,0) = Vexp 



z 



igj t a A a (y)dy 



n 



(2.5) 



is the Fock-Schwinger phase factor (termed in [23] the color "connector"), path-ordered 
along the straight line joining the points and z, to preserve gauge invariance. Note that 
the scale /Xq, called the normalization scale of the pion DA, is related to the ultraviolet 
(UV) regularization of the quark-field operators on the light cone whose product becomes 
singular for z 2 = 0. The derivation of hadron distribution amplitudes at finite momentum 
transfer is outside perturbative QCD; it requires nonperturbative approaches, like QCD sum 
rules with vacuum condensates [19], or lattice calculations. The latter method has recently 
reached a high level of precision in calculating the first and second moment of the pion DA 
[24, 25, 26], though the calculation of the fourth moment is still pending, while recently 
a prediction for its value was worked out [27] by combining the CLEO data [28] and the 
lattice calculation for the second moment. Here, an alternative method proposed by Braun 
and Miiller [29] may be proven useful. On the other hand, QCD sum rules incorporate 
information about the non-trivial structure of the QCD vacuum and, therefore, constitute a 
useful analytic tool for determining hadron distribution amplitudes using local constraints 
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on their moments (inverse-scattering method). Chernyak and Zhitnitsky [19] extracted the 
lowest few moments (£ 7V ) = (p n (x) (2x—l) N dx of <p n using correlators of two axial currents 
with local operators containing iV covariant derivatives [1, 5]. However, it was pointed out 
later [15, 30, 31] that QCD sum rules with local condensates strongly emphasize the endpoint 
region, because the condensate terms are strongly peaked just at these endpoints x = 0,1, 
where one of the quarks has a vanishing virtuality. Therefore, the authors of Refs. [15, 30, 31] 
suggested to introduce nonlocal condensates that can account for the possibility that vacuum 
quarks may have a finite virtuality (a nonzero average transverse momentum). As a result, 
Feynman-type configurations, in which one of the quark carries almost the entire available 
momentum while the spectators are "wee" (with almost vanishing virtualities), are separated 
out from the pion DA and treated separately in a non-factorizing (soft) contribution to the 
form factor (for a discussion of the Feynman mechanism in this context, see [32]). 

More recently, a QCD sum-rule analysis was carried out [17] which employs nonlocal 
condensates parameterized for simplicity in terms of a single mass-scale parameter A 2 = (kj_) 
with values in the range [33, 34, 35, 36] 



Xz = 0.35 - 0.55 GeV 2 



(2.6) 



This way, the first ten moments of ip n at the initial scale /Xq = 1-35 GeV 2 were determined 
with the following values (theoretical errors in parentheses) [17] 



(?) = 0.265(20), (?) = 0.115(12), (?) « 0, (N = 6, 8, 10) 

3/2, 



(2.7) 



Recasting the pion DA in terms of the Gegenbauer polynomials C n (2x — 1), which are the 
one-loop eigenfunctions of the ERBL kernel [2, 4], one finds 



^(x,// 2 ) = ip as (x) 



l + ^a 2n (/i 2 )C 2 3 / 2 (2x-l) 



n>l 



y^Q2n(/i S 



n>l 



where the asymptotic pion DA has the form 

</ s ( x ) = 6x(l - x) 



:2.8) 



(2.9) 



The values of the Gegenbauer coefficients, encoding the nonperturbative dynamics, of some 
characteristic pion DAs are given in Table I, while the corresponding profiles are displayed 
in Fig. 1. 

It is worth noting that the rate of convergence <p n (x, fi 2 ) — > y9 as (x), or a 2n (Q 2 ) —> 0, at 
Q 2 — > oo is not fast and is determined in the one-loop approximation by the logarithmic law 



G2n(/-t 2 



&2n(/^o) 



'2n 



(2.10) 



DA 


02(1 GeV 2 ) 


o 4 (l GeV 2 ) 


GeV 2 ) 


Shape 


Asy 


0.00 


0.00 


3.00 


convex 


CZ 


0.56 


0.00 


4.68 


double-humped, end-point enhanced 


BMS 


0.20 


-0.14 


3.18 


double-humped, end-point suppressed 



TABLE I: Main characteristics of three different DAs, abbreviated by acronyms: CZ stands for 
Chernyak-Zhitnitsky [19] and BMS for Bakulev-Mikhailov-Stefanis [17]. Here a 2 and 04 are the 
second and fourth Gegenbauer coefficients, respectively, whereas IZ% denotes the inverse moment 
of <p n . 
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FIG. 1: Comparison of selected pion DAs labeled by obvious acronyms: ip aB (dotted line), <pcz (dashed line) 
[19], and ^bms ( s °lid unc ) [17], defined by the coefficients 02 and 04 in Table I. All DAs are normalized at 
the same scale /z(j ~ 1 GeV 2 . 

Here 73° = 0.62, 7^ = 0.90, whereas all other anomalous dimensions 72° > 1 for n > 
3. Numerical estimates show that if one has at some typical initial scale fi 2 ~ 1 GeV 2 
a coefficient a 2 (l GeV 2 ) = 0.25, as indicated by the CLEO data [28] and recent lattice 
simulations [37, 38] , its value would become 3 times smaller only at the tremendous scale of 
fi 2 ~ 75400 GeV 2 , which is certainly outside the reach of any experimental measurement. 
The situation at intermediate momentum transfers 20 GeV 2 > Q 2 > 4 GeV 2 is more delicate. 
In this region, factorization partially fails owing to a quite sizable contribution of the soft 
part which is non-factorizable. Hence, one has to calculate this part of the pion form factor 
using either phenomenological models [7, 8, 39, 40], or by employing some nonperturbative 
concepts, like the method of QCD sum rules (SR) [13, 41, 42]. Still another option is to apply 
the local quark-hadron duality approach [14, 43, 44]. On the other hand, at the one-loop 
level and at asymptotically large Q 2 , the pion form factor turns out to be [45, 46] 

1 

F pert (g2) = STT^Xg , , )2 ^ 2 = t ^Q*) ^ 

9Cr J x 



The precise value of Q 2 at which this perturbative expression should start to prevail cannot 
be predicted (determined) accurately. The estimates for the cross-over momentum scale 
range from 100 GeV 2 [7, 8, 47] down to values around 20 GeV 2 [9, 10, 32]. But even this 
latter relatively small momentum is hopelessly far away from the capabilities of any operating 
or planned accelerator facility. 

III. QCD VACUUM AND THE PION FORM FACTOR 

In view of these facts, one may look for alternative methods to calculate the pion form 
factor in the experimentally accessible region of momentum transfer. Even after the com- 
missioned 12 GeV upgrade for Continuous Electron Beam Accelerator Facility (CEBAF), 
the expected high-precision experimental data will not be sufficient to enter the perturba- 
tive regime. Hence, instead of predicating to the elusive hope for still higher energies in the 
future, we discuss an alternative approach, based on three-point QCD SR to the pion form 
factor [41, 42]. One of the advantages of this technique is that the shape of the pion DA is 
irrelevant, reducing this way the theoretical uncertainty. 
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The construction of the Borel SR for the pion form factor from the three-point AAV 
correlator 

(3.1) 



JJd'xd%e^-^H0\T[j+(y)J^x)J 5a (0)]\0} 



where q corresponds to the photon momentum (q 2 = —Q 2 ) and p 2 is the outgoing pion 
momentum, was described in detail in [43, 48] for the case of local condensates and in [13] 
for the non-local condensate (NLC) case. In this correlator, J^(x) = e u u{x)^f^u{x) + 
edd(x) r y tJ 'd(x) is the electromagnetic current, while J^ a {x) = d(x)j5 r y a u(x) and J^(x) = 
w(x)757 'pd(x) are axial- vector currents, where e u = 2/3 and = —1/3 are the electric 
charges of the u and d quarks, respectively. Herewith one obtains the following SR: 



«0 SQ 



dsi ds 2 p 3 (si,S2,Q 2 )e 



-(S!+S 2 )/M 2 



+ $ G (Q 2 ,M 2 ) + <S> m (Q 2 ,M 2 ), (3.2) 







where the quark condensate contribution 



$ (9 - 9> (Q 2 , M 2 ) = $ 4Q (Q 2 , M 2 ) + $ 2V (Q 2 , M 2 ) + <$> qAq {Q\ M 2 ) 

contains the four-quark condensate (4Q), the bilocal vector-quark condensate (2V), and the 
antiquark-gluon-quark condensate term (q~Aq). The gluon-condensate contribution to SR 
(3.2) is presented by the Qq{Q 2 ,M 2 ) term. The graphical illustration of the corresponding 
diagrams is shown in Fig. 2. The perturbative three-point spectral density entering the SR 
above reads 



pf( Sl ,s 2 ,g 2 ) + ^4^ApW( Sl , S2 ,g 2 ) 



47T 



Recall that the leading-order spectral density 



pf\s 1 ,s 2 ,t) 



Air 2 



dt 2 



?d*_ 
3 dt 3 



1 



- yj(s! + s 2 + tf 



(3.3) 



(3.4) 



4sis 2 

has been calculated in the early eighties [43, 48], whereas the explicit (but too compli- 
cated to show it here) expression of the analogous next-to-leading order (NLO) version 
Apg (si, s 2 , Q 2 ) has been obtained quite recently in [49]. Note that the contribution from 
higher-resonances (HR), -Fhr, is usually modeled with the help of the same spectral density 



Phr(si, s 2 ) = [1 - 0(si < s )9(s 2 < s )} p 3 (si, s 2) Q 7 



(3.5) 



and using the continuum threshold parameter Sq. 

A novelty of the present investigation is the use of a running coupling that shows by 
construction analytic behavior in Q 2 , i.e., has no Landau singularities. This is done by 
adopting the arguments and techniques used in our previous works in [10, 44]. To be specific, 
we will use the one-loop analytic Shirkov-Solovtsov coupling [50] 



a s (Q 2 



Air 



A 2 



b Un(Q7A 2 QCD ) Q 2 -A 2 QCD 



(3.6) 



with bo = 9 and Aqcd = 300 MeV. The interested reader can find more about this subject 
in the recent reviews [51, 52, 53]. The nonperturbative terms $g an d &(qq) in the local- 
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(4Q) 



(2V) 



(qAq) 



(G) 



FIG. 2: Nonperturbative contributions [13, 43, 48] to the QCD SR, (cf. Eq. (3.2)). Here we show only 
typical for each subclass diagrams, whereas the complete set includes also mirror-conjugated diagrams (for 
the 4Q and qAq subclasses) and diagrams with permutations of the gluon-lines insertions (subclass G). 



condensate case are well-known [43, 48]: 

(a s GG) 



$£ C (M 2 



12ttM 2 



M 



104 A c 
M 4 



2Q 2 
13 M 2 



(3.7) 



The standard QCD SRs [43, 48] for the pion form factor are based on these local expressions 
which show a wrong scale behavior at large Q 2 . This means that the quark contribution 
contains both a linearly growing term as well as a constant one, while the gluon contribution 
is just a constant (see Table II). At the same time, the first, i.e., the perturbative, term on 
the right-hand side of Eq. (3.2) behaves at large Q 2 as s /Q 4 or M 2 /Q A . For this reason, 
the SR becomes unstable for Q 2 > 3 GeV 2 . Hence, in order to be able to extract predictions 
from a SR of this sort in the Q 2 region between 3 — 10 GeV 2 , we have first to improve the 
quality of the SR. 

TABLE II: <5 2 -behavior of the nonperturbative contribution in different QCD SR approaches. Here 
ci, C2, C3, C4 are dimensionless constants (not depending on Q 2 ). The abbreviations used are: LD 
for local duality, LO for leading order, and NLO for next-to-leading order, while A 2 and M 2 denote 
the vacuum quark nonlocality parameter and the Borel parameter, respectively. 



Approach 


Accuracy 


Condensates 


Q 2 -behavior of $ope 


Standard QCD SR [43, 48] 


LO 


local 


ci + Q 2 /M 2 


QCD SR with NLCs [13] 


LO 


local + nonlocal 


(c 2 + Q 2 /M 2 ) L-ca&fyM* + C4 ) 


LD QCD SR [41, 47, 49] 


NLO 


no 





Here 


NLO 


nonlocal 


(ci + Q 2 /M 2 ) e ~ c ^ 2x V MA 



The fact that the condensate contributions in the SR for the pion form factor are con- 
stant, or even growing with Q 2 , is somewhat surprising because exactly the corresponding 
diagrams should actually generate decreasing contributions as Q 2 increases. But recall that 
the diagrams generating the condensate contributions differ from the ordinary Feynman 
diagrams of QCD perturbation theory. They result from the replacement of some of the 
propagators by constant factors which represent condensates. For example, the quark prop- 
agator (T(q(z)q(0))) is substituted by the quark condensate (g(O)g(O)). As a result, instead 
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of obtaining a Q 2 - dependent contribution, one gets a constant one. The dependence on Q 2 
appears when one calculates the contributions of higher-dimension operators of the type 
(q(0)D 2 q(0)) , (q(0)(D 2 ) 2 q(0)) etc., entailed by the Taylor expansion of the original nonlocal 
condensate (NLC), where (q(0)q(z)) is the nonperturbative part of the quark propagator. 
The resulting total condensate contribution decreases for large Q 2 . However, each term of 
the standard OPE has the structure (Q 2 /M 2 ) n , and one should, therefore, resum them to 
get a meaningful result. Our strategy is to avoid the original Taylor expansion and deal 
instead directly with the NLC. This leads to a modified diagram technique involving new 
lines and vertices that correspond to the NLC. Then, the simplest contribution to &{q q ) is 
due to the vector condensate M M (cf. (A.2)); viz., 



A$ y (M 2 ,Q 2 



8A 
M 4 



2 + 



Q 2 



M 2 



\2 



exp 



Q 2 K 



2M 2 (M 2 - A 2 ) 



(3. 



Obviously, this term indeed vanishes for large Q 2 , as expected. Moreover, the larger the 
nonlocality parameter A 2 , the faster this contribution decreases with Q 2 . The value Q 2 at 
which this decrease starts, strongly depends on the value of M 2 . Adopting for the nonlocality 
A 2 = 0.4 GeV 2 , one finds for M 2 = 1, 1.5, 2 GeV 2 , Q 2 = 2, 6, 13 GeV 2 , respectively. 




q[x) 



FIG. 3: Illustration of the distance dependence of the three-point condensate in terms of the mutual 
inter-parton separations. 



A first attempt [13] to generalize the QCD SR to the NLC case employed the specific model 
(3.9) for the three-point quark-gluon-quark NLC (with the NLC Mi(x 2 , y 2 , (x — y) 2 ) given by 
Eq. (A.2)) and assuming a nonlocality only with respect to two inter-parton separations, say, 
x 2 and (x — y) 2 , out of the three possible separations x 2 , y 2 , and (x — y) 2 (see Fig. 3). Hence, 
this approach is only partially nonlocal. The following parametric functions (A = A 2 /2) were 
used: 

ff R (a, 0,7) = 8 (a- x a A) 8 (J3 - x i2 A) 8 ( 7 - x i3 A) , (3.9) 
/OA 0.4 \ 
Xij = 1 0.4 . 
\ 0.4 0.4/ 

The zero elements in the matrix Xij indicate the absence of nonlocality effects either for 
the quark- ant iquark separation y 2 (i = l,j = 2) or for the antiquark-gluon separation x 2 
(i = 2, 3 and j = 1) — see Fig. 3. Therefore, also this approach partially suffers from the 
same shortcomings as the standard QCD SR. 

Fortunately, the NLC contributions to the pion form factor taken into account in the SR 
of [13], have a model-independent form, allowing us to use them in connection with more 
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improved versions of the quark-gluon NLC To be more specific, we apply here QCD SR 
employing the minimal (A. 3) [15, 16, 17] and the improved (A. 4) [18] Gaussian models of 
NLC. The corresponding parameters of the NLC are listed in Appendix A. An advantage of 
employing NLCs is that their use considerably enlarges the region of applicability of the QCD 
SR to momenta as high as 10 GeV 2 . In addition, we extend the accuracy of the perturbative 
spectral density to the NLO level by including contributions of 0(a s ) [49] — in contrast to 
previous works [13, 43, 48] in which only a LO perturbative spectral density was taken into 
account. It turns out that the NLO contribution influences the prediction for the pion form 
factor, calculated with the described method, reaching the level of 20%. 

IV. PION FORM FACTOR ANALYSIS 
A. Calculation and predictions 

The strategy to further process the obtained SR is standard: At each fixed value of Q 2 , SR 
(3.2) gives us the pion form factor F n (Q 2 , M 2 , sq) as a function of two additional parameters, 
notably, the auxiliary Borel parameter M 2 and the continuum threshold Sq. The parameter 
s can be interpreted as the boundary between the pion and the higher resonances (A±, ir', 
etc.). We assume that Sq should not be lower than the middle point, 0.6 GeV 2 , of the interval 
between m 2 = and m 2 Ai ps 1.6 GeV 2 . 

The specific value of s (Q 2 ) at a given value of Q 2 is determined by the condition implied 
by the minimal sensitivity of the function Ft T (M 2 ,sq) on the auxiliary parameter M 2 in 
the fiducial interval of the SR. We derive these intervals and the values of the pion decay 
constant f n from the corresponding two-point NLC QCD SR, employing both the minimal 
and the improved model for the vacuum nonlocality — see Table III for details. Note here 
that the value of the Borel parameter M 2 in the three-point SR roughly corresponds to the 
Borel parameter in the two-point SR, having, however, twice its magnitude: M t 2 hree point = 
2-^two-point) as gi ven hi Table III. In the left panel of Fig. 4, we show how the scaled pion form 
factor Q 2 F^Q 2 , M 2 , sq) depends on M 2 for three different values of the threshold parameter: 
s = 0.65, 0.75, and 0.85 GeV 2 at Q 2 = 5 GeV 2 . As a rule, the higher the value of sq, the 
larger the form factor because the perturbative input increases. 

Using the root-mean-square deviation x 2 {Q 2 i s o), given by Eq. (B.l), we determine that 
continuum threshold Sq R (<5 2 ) which minimizes the dependence of the right-hand side of (3.2) 
on the Borel parameter M 2 e [M 2 , M 2 ] at each value of Q 2 . For an illustration, we refer to 
the right panel of Fig. 4. 

As one sees from the right panel of Fig. 4 the long-dashed line — which corresponds to 
the minimal NLC model at Q 2 = 5 GeV 2 — has no minimum in the relevant so interval. 
Therefore, we cannot extract a reliable continuum-threshold value Sq R (Q 2 > 4 GeV 2 ) for 
this model using the root-mean-square deviation criterion. 



TABLE III: Boundaries of the fiducial intervals M|/2 with M 2 £ [M 2 /2, M 2 /2] and values of the 
pion-decay constant (/„-) used in the two-point sum rules with nonlocal condensates in connection 
with the minimal and the improved Gaussian model for the nonlocality. 



Model 


u 


Ml 


M\ 


Minimal [17] 


0.137 GeV 2 


1 GeV 2 


1.7 GeV 2 


Improved [18] 


0.142 GeV 2 


1 GeV 2 


1.9 GeV 2 



9 



. 65 

. 6 
0.55 

. 5 
.45 

0.4 
0.35 

. 3 



Q 2 F„(Q 2 ,M 2 ,s ) [GeV 




FIG. 4: Left panel: Dependence of the pion form factor Q 2 F 7T (Q 2 ) at Q 2 = 5 GeV 2 on the auxiliary Borel 
parameter M 2 in the improved NLC model. The solid blue line corresponds to sq = 0.75 GeV 2 , whereas 
the dashed lines refer to s = 0.65 GeV 2 (upper curve) and so = 0.85 GeV 2 (lower curve). Right panel: 
Root-mean-square deviation x 2 (Q 2 i s o) f° r the minimal NLC model at Q 2 = 3 GeV 2 (long-dashed blue line) 
and at Q 2 = 5 GeV 2 (short-dashed red line) . The improved model of NLC is shown as a solid blue line at 
Q 2 = 5 GeV 2 . 

Note, however, that the values min[x 2 (Q 2 , s)] and x 2 {Q 2 , s o = Sq D '^\Q 2 ) ^ 0.63 GeV 2 ) 

s 

are very close to each other, the relative difference for Q 2 = 4 — 10 GeV 2 being of the order 
of 10 — 15%*. For this reason, we will use in the minimal model Sq R (Q 2 ) = Sq D (Q 2 ) 
as the continuum threshold. In contrast, in the improved case (the solid line in the right 
panel of Fig. 4), the root-mean-square deviation has a minimum approximately at the same 
point Sq R (Q 2 ) ~ 0.75 GeV 2 for any Q 2 value within the considered interval. The continuum 
thresholds Sq R (C} 2 ) for the minimal (dashed line) and the improved (solid line) NLC model, 
obtained this way, are shown in Fig. 5. 

On that basis, we can define the SR result for the pion form factor as the average value 



i 

0.8 
0.6 
. 4 
0.2 



s (Q 2 ) [GeV 2 



Q 2 [GeV 21 



10 



FIG. 5: Continuum threshold sq(Q 2 ) [GeV 2 ] for the minimal (dashed line) and for the improved (solid line) 
NLC model. 



Here Sq D '^ (Q 2 ) is the standard Local-Duality prescription for the continuum threshold, see the discussion 
after Eq. (4.5). 
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of the right-hand side of (3.2) with respect to the Borel parameter M 2 £ [M 2 , M?]: 

F™(Q 2 ) = -^TZW C F(g2 ' S ° SR(g2)) ^ • 

Our main predictions for the minimal and the improved NLC model, using in both cases 
A 2 = 0.4 GeV 2 , are shown in Fig. 6 in the form of two bands, each one corresponding to the 
particular NLC model used, with the width of each band denoting the inherent theoretical 
uncertainties of the underlying QCD SR method. The band within the dashed lines contains 
the predictions for the minimal model. Its counterpart for the improved model is limited by 
the solid lines. For the central curves (dashed — minimal model; solid — improved model) in 
Fig. 6, as well as in both panels of Fig. 7, we used the following interpolation formulas: 

F™ in (Q 2 = x GeV 2 ) = e- 1 ' 402 * ' 525 (l + 0.182* + °~f^f) , (4.2a) 
*SLp(Q a = x GeV2 ) = e- 1 - 171 * ' 536 (l + 0.0306 x + °'f+f ) , (4.2b) 

valid for Q 2 £ [1,10] GeV 2 , i. e., for x £ [1,10]. The two broken vertical lines in Fig. 6 
denote the strict fidelity window of the NLC QCD sum rules. We have limited this window 
from above by the requirement that the predictions obtained with both NLC models have 
an overlap. As one sees, the central curve for the minimal model starts departing from the 
band of the improved model just around 7 GeV 2 . But the predictions remain useful (though 
less accurate) even at higher Q 2 values close to 10 GeV 2 . Our results compare favorably with 
the lattice calculation of [56], shown as a monopole fit with associated error bars between 
the two thick lines at lower Q 2 . We observe a similarly good agreement with the existing 
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FIG. 6: Scaled pion form factor Q 2 F 7T (Q 2 ) for the minimal (dashed blue lines) and for the improved 
(solid blue lines) NLC model using X 2 . = 0.4 GeV 2 in comparison with experimental data of the Cornell 
[11, 54, 55] (triangles) and the JLab Collaboration [12] (diamonds). The shaded band delimited by the 
dashed lines corresponds to the minimal model, whereas the analogous band bounded by solid lines represents 
the uncertainty range of the improved model. The recent lattice result of [56] is shown as a monopole fit 
with error bars between the two thick lines at lower Q 2 . The two broken vertical lines denote the region, 
where the influence of the particular Gaussian model used to parameterize the QCD vacuum structure in 
the NLC QCD sum rules is not so strong. 
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experimental data of the Cornell [11, 54, 55] (triangles) and the JLab Collaboration [12] 
(diamonds). 

Figure 7 shows our results in comparison with the predictions obtained with various 
theoretical models (left panel), whereas the right panel compares our results with predictions 
from AdS / QCD models. Specifically, in the left panel of Fig. 7 we display as a short-dashed 
green line the prediction from [57], whereas the long-dashed green line is from [58]. The short 
heavy solid red line at low Q 2 represents the standard QCD SR result with local condensates 
[43, 48], while the dash-dotted red line gives the more recent estimate of the Local Duality 
QCD SR of [47]. Note here that the results of [10], which are not shown in this figure, are 
approximately 20% higher than those represented by the dash-dotted red curve. This is due 
to (i) the inclusion of the 0(a 2 ) correction (10%) and (ii) the uncertainty of the matching 
procedure, used in [10], (10%) (see the discussion in Sect. IV B and the graphics in Fig. 9). 
Finally, the dotted line below 4 GeV 2 displays the result extracted from the model of [59] 
which uses the Bethe-Salpeter equation. 

In the right panel of Fig. 7 we collect recent estimates from the following AdS/QCD 
models: The heavy long-dashed line displays the result obtained with an AdS/QCD pion 
DA [60]. The dash-dot-dotted green line represents a Hirn-Sanz-type holographic model, 
discussed in [61]. The dotted green line provides the prediction of the AdS/QCD soft model 
[62]. Finally, the two top broken red lines show the results obtained from the improved 
soft-wall (heavy dash-dotted line) and the hard- wall (fine dash-dotted line) AdS/QCD back- 
ground approach [63], respectively. One sees that the dotted and the dash-dot-dotted green 




2468 10 2468 10 



FIG. 7: Comparison of theoretical predictions for the scaled pion form factor Q 2 F 7T (Q 2 ) obtained with 
a variety of theoretical methods and models. Our estimates are shown in both panels as shaded bands 
corresponding to the minimal NLC model (dashed blue lines) and the improved NLC model (solid blue 
lines). In both cases X 2 = 0.4 GeV 2 has been used and the associated error range is indicated by the width 
of the shaded bands. Left panel: The first short-dashed green line above to the dash-dotted red one is from 
[57] and is based on a large- N c Regge model. The long-dashed green line is from [58] , obtained by including 
radiative and higher-twist effects within the framework of resummcd pQCD. The short solid red line at low 
Q 2 shows the result from the standard QCD SR with local condensates [43, 48]. The dash-dotted red line 
denotes the estimate from [47] which was derived from Local Duality QCD SR. The dotted line represents 
the model of [59], based on the Bethe-Salpeter equation. Right panel: The heavy long-dashed line shows the 
result obtained with an AdS/QCD pion DA in [60]. The dash-dot-dotted green line shows a Hirn-Sanz-type 
holographic model result, discussed in [61]. The dotted green line gives the prediction of the AdS/QCD 
soft-wall model [62]. The heavy and the fine dash-dotted red lines show, respectively, the predictions derived 
from the improved soft and hard wall AdS/QCD background approach [63]. 
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FIG. 8: Contributions to the pion form factor from the O (1) and O (a s ) spectral-density terms (LO and 
NLO). The specific contributions are denoted by labels referring to their particular origin: 4-quarks (scalar) 
condensate — (4Q), bilocal vector quark condensate — (2V), tri-local quark-gluon-quark condensate — (qAq), 
and gluon condensate — (G) in the minimal (left panel) and the improved (right panel) model. 



lines are in compliance with our predictions, favoring the results of [61], as these leave some 
space to add radiative corrections. 

Figure 8 serves to illustrate the origin of the differences between the minimal (dashed 
line — left panel) and the improved (solid line — right panel) Gaussian NLC model with respect 
to their corresponding predictions for the pion form factor. The difference between them 
can be traced back to the quark-gluon-quark contribution $>qAq(Q 2 , M 2 ) which is different 
for each model (top lines in both panels below the final results). This, in turn, has influence 
on the continuum threshold Sq(Q 2 ) entailing changes in the corresponding LO and NLO 
terms, albeit mild ones. This figure exhibits in detail how the various contributions to the 
pion form factor are accumulated. Each curve from bottom to top is the sum of all previous 
terms, while the heavy lines on the top denote the full result for each model. 

Note that the calculations were made for a nonlocality parameter A 2 = 0.4 GeV 2 , a value 
receiving support from a recent comprehensive analysis [64, 65] of the CLEO data on the 
pion-photon transition. Using higher values of this parameter, would entail a decrease of 
the pion form factor owing to a stronger influence of the nonlocality effects. Evidently, 
using smaller values of A 2 would cause an increase of the pion form factor. The predictions 
shown in our figures, discussed above, provide further support for a value of A 2 in the range 
compatible with the result extracted from the CLEO-data analysis. 



B. Local duality approach 

Sum rules based on local duality have no condensate contributions due to the M 2 — > 00 
limit. On the other hand, the determination of the threshold sq is not possible using this 
method, in contrast to QCD SR. One could, nevertheless, try to extract it from (local duality) 
sum rules for the pion-decay constant f n , but this would correspond to a value of the pion 
form factor at a small Q 2 value because of the Ward identity [47]. 

As originally argued in Refs. [7, 41, 43, 48, 66], the dominant contribution to the pion 
form factor at low up to moderate values of the momentum transfer Q 2 < 10 GeV 2 originates 
mainly from the soft part that involves no hard-gluon exchanges but may be attributed to 
the Feynman mechanism. Using the Local Duality (LD) approach to calculate the soft 
contribution, it is assumed that the pion form factor is dual to the free quark spectral 
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density [14, 43]. Then, staying within the (I + l)-loop order^, one has 



F LD;( 0(g2) = F LD ; (0 (g2; (4 3a) 

F^ l \Q\S) = 1 [ S [ S p!\s 1 ,s 2 ,Q 2 )ds 1 ds 2 , (4.3b) 

where Sq ''''^^ 2 ) is the LD effective threshold parameter for the higher states in the axial 
channel and the three-point (I + l)-loop spectral density is p$ (si, S2, Q 2 )- In leading order, 
we know pf\si, S2,Q 2 ) from Eq. (3.4), so that 



F ^ s) = j_^^ ] j_^_y „ 4) 

The LD prescription for the corresponding correlator [14, 67] implies the relations 

So "«>) =4^/1 and C i( "(0) , (4-5) 

1 + a s {Q )/7T 

where Ql is of the order of Sq '^^). This prescription is a strict consequence of the Ward 
identity for the AAV correlator due to the vector-current conservation. In principle, the 
Q 2 dependence of the LD parameter Sq D (Q 2 ) (4.3) should be determined from the QCD 
SR at Q 2 > 1 GeV 2 . As we have already explained in Sec. Ill, the standard QCD SR 
becomes unstable at Q 2 > 3 GeV 2 because of the appearance of terms in the condensate 
contributions linearly growing with Q 2 [41, 42]. For this reason, this dependence was known 
only for Q 2 < 3 GeV 2 and, therefore, most authors usually used the constant approximation 
Sq (Q 2 ) — s o D ( ' ^(0)5 like in [10, 43, 44, 49], or a slightly Q 2 -dependent approximation 

LD;(1) 



'o 



(Q 2 ) ^4n 2 f 2 /(l + a s (Q 2 )/*), like in [47]. 



Lacking profound knowledge about the exact structure of the NLO spectral density 
p£\si, s 2 , Q 2 ), two of us (A.B. and N.G.S.) with collaborators have suggested in [10] to 
use for the full pion form factor the information about the factorizable part, F^ CiCI:> ' < ' 2 \Q 2 ), 
which is computable within perturbative QCD (pQCD) (here evaluated at the two-loop 
order). But it should be noted that the pQCD term has the wrong limit at Q 2 = 0, calling 
for a correction of this behavior in order to maintain the Ward identity (WI) ^(0) = 1. To 
achieve this goal, we apply a matching procedure, introduced in [10], namely, 

F WI;(2) (g2) = F LD,(0) (g2) + ^Q 2 ^ _ \ F pQCD,( 2)(g2) (4 g } 

with ~ 0.6 GeV 2 . This approximation was used to 'glue' together the LD model for 
the soft part, Fn D ^°\Q 2 ) (which is dominant at small Q 2 < 1 GeV 2 ), with the perturbative 
hard-rescattering part, F^^'^iQ 2 ) (which provides the leading perturbative 0(a s )+0(a 2 ) 
corrections and is dominant at large Q 2 3> 1 GeV 2 ), in such a way as to ensure the validity 
of the Ward identity Fj^'^O) = 1. In order to test the quality of the matching prescription 



t The triangle diagram at the one-loop level docs not include any radiative corrections from the very begin- 
ning. For this reason, the 0(aJ,)-correction for this diagram appears first at the (I + l)-loop order. 
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given by Eq. (4.6), we propose to compare it with the LD model (4.3) evaluated at the 
one-loop order (i.e., in the 0(a s )-approximation [47, 49]). To this end, we construct the 
analogous 0(a s )-model 



F WI;(l)/ ( o2N _ F LD,(0)(Q2\ i a s{Q 2 ) 2 Q 2 Sq D ' (0) (0) 

where we made use of the asymptotic form factor [45, 46] 

FP QCD,(i )rr - _ Stt/^Q 2 ) _ 2s^°\o) 
W)~ Q2 ~ n Q2 , 

implying the same prescription for the effective LD threshold as in [47] : 

W = T^Sfe' <4 ' 7C) 

It is worth noting that the model F^ 1 ' (Q 2 ), following from the matching procedure (4.7a) 
suggested in [10], works quite well, albeit it was proposed without the knowledge of the 
exact two-loop spectral density, which became available only somewhat later [49]. Recall 
that the key feature of the matching prescription is that it uses information on F n (Q 2 ) in 
two asymptotic regions: 

1. Q 2 0, where the Ward identity dictates F n (0) = 1 and hence F n (Q 2 ) ~ i^ D ' (0) (Q 2 ), 

2. Q 2 -> oo, where F^Q 2 ) ~ F. pQCD ' (1) (Q 2 ) 

in order to join properly the hard tail of the pion form factor with its soft part. Numerical 
analysis of (4.7a) shows that the applied prescription yields a pretty accurate result, with 
a relative error varying in the range 5% at Q 2 = 1 GeV 2 to 9% at Q 2 = 3 — 30 GeV 2 . 
The graphical comparison of F™'' {1) (Q 2 ) (solid blue line) with the LD result (4.3) (black 
dots), employing the two-loop spectral density [49], is displayed in the left panel of Fig. 9. 
This figure also shows the purely perturbative part Q 2 F^ CD '^\Q 2 ) (dashed red line) — not 
corrected by the factor [Q 2 /{2sq + Q 2 )} 2 to the right low-Q 2 behavior. * [In similar context 
a matching procedure, as that described above, was also used in a related work [58] on the 
electromagnetic pion and kaon form factor including radiative and higher-twist effects.] 

Once the spectral density pg (s\, s 2 , Q 2 ) was calculated [49], it made it possible to improve 
the representation of the LD part in (4.7a) by taking into account the leading 0(a s ) cor- 
rection in the electromagnetic vertex. On that basis, we suggest the following improved WI 
model for F^ (Q 2 , s^iQ 2 )): 

F^iQ 2 ,S) = F^\Q 2 ,S) 

^(asms^s^ 2 QCD , ( 1)(g2) (_£_y| (48) 

+ 4vr 2 / 2 \ tt \2S + Q 2 ) + ^ W) \2S + Q 2 ) J ' ^ 

We explicitly display the dependence on the threshold S in Eq. (4.8) — the aim being to 
apply it later on with S = Sq D ~ cR (Q 2 ), the latter value being extracted by comparing the 



$ For this reason, this expression tends to the finite value 0.21 GeV 2 at Q 2 — > and does not vanish. 
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FIG. 9: Comparison of the exact NLO LD result for the pion form factor F^ (Q 2 ) (black dots) 



with the models F™ ;(1 \q 2 ) ((4.7a), left panel) and F^> (Q 2 , s^ 1 (Q 2 )) ((4.8), right panel) : 

displayed in both panels as solid blue lines. The asymptotic perturbative prediction F^ CD '^ 1 \q 2 ) 
is also shown dashed red line. 



WI ;(l)//->2 LD;(1) /^- ) 2 N i 



NLC QCD SR results with the LD approximation. Note that this model has the right limit 



at Q 2 — > 0, provided one sets S = So A ) (Q 2 ) = 4vr 2 / 2 /(l + a s (Q 2 )/n). Indeed, one has 
i^ LD;(0) (0,S) = S/(4ir 2 f 2 ) by virtue of the fact that the term F% QCD ' {1) (Q 2 ) is canceled by 
the factor [Q 2 /(2S + Q 2 )] 2 in this limit, so that the net result is § 



pWI;(l) 
7r;imp 



(0. 



,LD;(1) 



(0)) 



,LD;(i) 

'0 



(0) 



47r2/ 2 



1 + 



7T 



The graphical evaluation of this new WI model in comparison with the exact LD result in 
the one-loop approximation is displayed in the right panel of Fig. 9. We can see from this 
graphics that the quality of the matching condition, given by (4.7a), is improved. Indeed, the 
relative error is reduced, varying between 4% at Q 2 = 1 — 10 GeV 2 and 3% at Q 2 = 30 GeV 2 . 
Proceeding along similar lines of reasoning, we construct the two-loop WI model 



^ WI;(2) (Q 2 ,^ D;(2) (Q 2 ))for the pion form factor to obtain 



F™* 2 \Q 2 ,S) 



+ 



F^°\Q 2 ,S) 

S fa s (Q 2 ' 
4vr 2 / 2 



7T 



2S 



2S + Q< 



F FAPT,(2) (g2 



Q- 



2S + Q< 



(4.9) 



where F^ APT '^ 2 \Q 2 ) is the analyticized expression generated from F^ CD '^ 2 \Q 2 ) using Frac- 
tional Analytic Perturbation Theory (FAPT) (see Refs. [68, 69, 70, 71, 72] and, in particular, 
[52, 73]) to get a result which appears to be very close to the outcome of the default scale set- 
ting (/iR = fip = Q 2 ), investigated in detail in [10]. The explicit expression forFj APT < (2) (Q 2 ) 
is given in Appendix C. This model provides the possibility to implement the perturbative 
QCD 0(a 2 )-resuIts for the pion form factor without performing an explicit three-loop cal- 
culation of the three-point spectral density, a welcome advantage, as this calculation is very 



§ One usually applies some "freezing" assumption Sq D,< ' 1 ' > (Q 2 ) = Ait 2 f 2 /(1 + a s (so)/7r) for Q 2 < sq — 
0.6 GeV 2 . In this case, the same "freezing" should be used in Eq. (4.8) for the argument of a s as well: 
&s{Q 2 ) —* Q-s(so) for Q 2 < so- Note that in order to treat the electromagnetic radius of the pion (i.e., the 
derivative of the pion form factor in the low Q 2 domain) correctly, one has to apply a different form of 
the OPE and the corresponding SR — sec for details in [41]. 



16 



1 


- sJO 2 ) \GeV 2 ] 




0.8 






. 6 






0.4 






0.2 




g 2 [GeV 2 ] 



2 4 6 8 10 



FIG. 10: Effective continuum thresholds s^ n °^(Q 2 ) (solid blue line) and s^mf n ff (<2 2 ) (dashed blue line) that 
approximate the NLC QCD SR results using the LD 0(a s (<5 2 ))-formulas. 

tedious. Moreover, the case of the one-loop approximation with the WI model (4.8) indicates 
that the relative error of this procedure is of the order of 10%. This level of accuracy allows 
us to estimate the weight of the 0(a 2 )-correction. We find that its relative contribution 
to the pion form factor is of the order of 10%, as has been estimated in [10, 52]. Hence, 
the relative error of our estimate is of the order of 1% — provided we take into account the 
0(a s )-correction exactly via the specific choice of so(Q 2 ), as done in (4.10). In order to 
use this formula, we only need to improve our knowledge about the effective LD thresholds 

^ D;(1) (Q 2 )- 

As we also stated at the beginning of this section, the problem of adjusting the continuum 
threshold parameter Sq D (Q 2 ) is of high importance for the LD approach. From our point 
of view, it should be determined by comparing the LD results with those derived via the 
Borel QCD SRs. In the previous section, we processed our NLC QCD SR and obtained the 
interpolation expressions (4.2a), (4.2b) which are valid for Q 2 e [1,10] GeV 2 . Now we can 
determine the corresponding effective thresholds s^hf (Q 2 ) an d Sq (Q 2 ) in the minimal 
and the improved NLC model, respectively, to find 

F%£> = ^ SR min (Q 2 ), F^(Q 2 ,s^(Q 2 )) = F^ p (Q 2 ). (4.10) 

The solutions of these equations, namely, Sq^^(Q 2 ) and s^^f (Q 2 ), are shown in Fig. 10; 
they can be represented in this range of Q 2 by the following interpolation formulas: 

SoMniQ 2 = x GeV 2 ) = 0.57 + 0.307 tanh(0.165x) - 0.0323 tanh(775x) , (4.11a) 
Sofmf (Q 2 = x Gey2 ) = °- 57 + °- 461 tanh(0.0954x) . (4.11b) 

We see that both thresholds turn out to increase only moderately with Q 2 . 

The results obtained for the pion form factor with our two-loop model, i.e., Eq. (4.9), 
and using the effective LD thresholds Sq^^(Q 2 ) and Sg^f (Q 2 ), are displayed in Fig. 11. 
We see from this figure that the main effect of the NNLO correction peaks at Q 2 > 4 GeV 2 , 
reaching the level of 3 — 10%, and can be estimated by taking recourse to the results of the 
FAPT analysis in [10, 52]. 

C. Estimation of uncertainties 

Even with our methodologically improved approach, there are still intrinsic uncertainties 
that influence our results. Therefore we consider in this subsection two potential sources 
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FIG. 11: We show as a narrow dash-dotted red strip the predictions for the pion form factor, obtained in 
the two-loop WI model, Eq. (4.9), using the minimal (left panel) and the improved (right panel) Gaussian 
model. The width of the strip is due to the variation of the Gegenbauer coefficients a,2 and 04 (needed 



to calculate the collincar part F% 



pQCD,(2) 



(Q 2 )) in the corresponding shaded bands for the pion DA in both 



models (indicated by the central blue lines). The normalization of the strip is also affected by the effective 
continuum thresholds So mS?(Q 2 ) panel) and 8q hwfCQ 2 ) (right panel). 



of uncertainties: (i) the choice of Aqcd and (ii) those uncertainties originating from our 
nonperturbative assumptions. 

• We performed our numerical estimates using a value of Aqcd = 300 MeV for 3 flavors 
of active quarks. Let us here discuss how our results change as we vary the value of 
this parameter. First, we note that the relative contribution of the NLO correction to 
the value of the pion form factor is of the order of 20%, as we already mentioned at the 
end of Section III. Next, the relative difference between a s (Q 2 , Aqcd = 300 MeV) and 
a s {Q 2 , Aqcd = 200 MeV) is of the order of 15% for Q 2 = 1 - 10 GeV 2 . Hence, varying 
the QCD scale parameter Aqcd i n the range 200 — 300 MeV induces uncertainties of 
our predictions of the order of about 3%. 

• As it was shown in [18], the minimal Gaussian model for the nonlocal quark condensate 
does not satisfy the QCD equations of motion and the transversality condition for the 
two-point correlator of the vector currents. For this reason, using the minimal model 
can induce some artificial term in the OPE part of the QCD SR. Bearing this in mind, 
the improved model was constructed in such a way as to minimize this unphysical con- 
tribution for the first five lowest-order conformal moments of the pion DA. Speaking in 
terms of a Taylor expansion of the nonlocal condensates, this improvement pertains to 
the correction of the low-dimensional local condensates — which represent coefficients in 
this expansion. The dark side of this improvement is that it could render the definition 
of higher-dimensional local condensates ill-defined relative to the minimal model. This 
is the reason why we refrained from applying the improved model to the extraction of 
the higher conformal moments of the pion DA and the high-momentum behavior of 
the pion form factor. As we can see from Fig. 8, in the case of the improved model the 
quark-gluon term (qAq) grows with Q 2 in the considered momentum region — in con- 
trast to the minimal model. We conclude from this that the qAq-term in the improved 
model is not well-defined in the high Q 2 region. One can realize from Fig. 7 that 
the difference between the minimal and the improved model tends to become larger 
relative to the accuracy of the SR at Q 2 ~ 7 GeV 2 . This means that the validity of 
the nonperturbative model used is lost at this point. Therefore, we should anticipate a 
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reduced accuracy of the obtained results on the pion form factor in the high-momentum 
regime. To overcome this problem, one should eventually set stricter constraints on 
the model using experimental data — a promising option after the planned upgrade of 
the CEBAF accelerator. Still another possibility comes from theory, notably, current 
conservation and the equations of motion. Such attempts are currently under scrutiny. 

V. DISCUSSION AND CONCLUSIONS 

In this work we have calculated the pion's electromagnetic form factor using a three-point 
QCD sum rule with nonlocal condensates. The main scope of our investigation was on the 
methodological side, though the obtained predictions compare well with both the existing 
experimental data and lattice simulations. The applied method offers the following key 
advantages: 

(a) The profile of the pion distribution amplitude is irrelevant. This removes a serious 
source of theoretical uncertainty entailed at the normalization point of the order of 
1 GeV and simplifies the computation because the technically demanding task of the 
NLO evolution of the pion DA to higher Q 2 values to compare with data is absent. 

(b) Because this type of approach tends to become unstable above Q 2 > 3 GeV 2 , we had 
to improve the quality of the sum rules considerably. This was achieved by dealing 
not with a Taylor expansion but directly with the nonlocal condensates, avoiding any 
local contribution. As a result, we were able to enlarge the region of applicability of 
the QCD SR towards momenta as high as 10 GeV 2 . 

(c) In order to minimize the transversality violations of the two-point correlator of vector 
currents, we considered in our analysis not only the minimal Gaussian model for the 
distribution of the quark average momentum in the vacuum, but also an improved 
version of it introduced earlier in [18]. 

(d) A spectral density was used in the dispersion integral of the sum rules which includes 
terms of 0(a s ). The influence of this next-to-leading-order (NLO) contribution to the 
pion form factor turns out to be quite important. 

(e) An analytic running coupling was used in the calculations pertaining to perturbation 
theory that helps surmount the problem of a Landau singularity at Q 2 = Aq CD . 

The phenomenological findings of our analysis and their implications are the following. 

(i) Our main predictions for F n (Q 2 ) were shown in Fig. 6 in comparison with the exist- 
ing experimental data from the early Cornell [11, 54, 55] and the more recent JLab 
Collaboration [12]. We found that the 0(a s )-contribution to the spectral density influ- 
ences the pion form factor at the level of 20%. This estimate contrasts with the result 
obtained in [47, 49], which was found to be somewhat larger. 

(ii) Inspection of Fig. 6 reveals that the central-line curve of the improved model lies 
within the error range of the minimal model up to values Q 2 ~ 7 GeV 2 , indicating 
a comparable quality of both models in this momentum region. We pointed out that 
the main differences between these two models for the vacuum nonlocality can be 
attributed to the different contributions related to the quark-gluon-quark contribution 
^WQ 2 ,M 2 ) (see Fig. 8). 
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(iii) In Fig. 7 we compared our predictions with a variety of theoretical models for the elec- 
tromagnetic pion form factor — including also results derived from holographic models 
based on the AdS/QCD correspondence. We found that the models of [61] and [62] 
are within the error bands of the minimal and the improved model, respectively. 

(iv) Applying a matching method of the factorized pion form factor, calculable within 
perturbative QCD, to its value at Q 2 = 0, subject to the Ward identity, we were 
able to model the pion form factor at the two-loop level without having recourse to 
the spectral density — appealing only to the local duality concept. This method was 
previously developed and applied successfully in [10] and was used quite recently with 
an appropriate modification in [58] in the calculation of the pion and the kaon form 
factors at the twist-three level. The key element in our approach is the extraction of 
the effective continuum thresholds Sq D (Q 2 ) from the results of the NLC QCD SR. 

(v) The results obtained with the local duality model, just described, were found to be 
inside the error band of the QCD sum rules with the nonlocal condensates (Fig. 11) 
providing support for the consistency of both approaches. Both methods yield predic- 
tions comparing well with the experimental data. 

(vi) Comparison of our predictions with those found with the LD approach [47] reveals that 
they are systematically higher than these. The reason for this difference is that the 
effective LD threshold Sq D (Q 2 ) has a well-defined value only in the small-Q 2 region. 
For higher Q 2 values, it is not firmly fixed. The authors of [47] proposed therefore a 
logarithmically increasing threshold 

47T 2 f 2 



l + a s (Q 2 )/7r' 

which is about 0.67 GeV 2 at Q 2 « 10 GeV 2 . In order to imitate the NLC QCD SR 
results within the LD approach, one needs to use Sq D (Q 2 = 10 GeV 2 ) = 0.87 GeV 2 . 
This means that the Sq D uncertainty in the region of Q 2 = 10 GeV 2 is of the order of 
20%. 

Bottom line: The use of three-point QCD sum rules with nonlocal condensates provides 
a reliable alternative to calculate the spacelike pion's electromagnetic form factor in that 
momentum region which is accessible to current and planned experiments. 
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APPENDIX A: PARAMETRIZATION OF THE NONLOCAL CONDENSATES 

We use, as usual in the QCD sum-rules approach, the fixed-point (Fock-Schwinger) gauge 
x^A^(x) = 0. For this reason, all connectors C(x,0) = V exp[—ig s J^ t a A c l 1 (y)dy fl ~\ = 1, 
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assuming for the integration contour a straight line going from to x. For the scalar and 
the vector condensates, we apply the same minimal model, as used in [16, 17]: 

(?(0)g(z)) = (qq) e -"*W 8 ; (?(0) 7/ ,?(*)> = A e"'^/ 8 , (A.l) 

where Aq = 2a s tt (qq) 2 / 81. The nonlocality parameter A 2 = (A; 2 ) characterizes the average 
momentum of quarks in the QCD vacuum and has been estimated in QCD SR [34, 74] 
and on the lattice [35, 36] to have a value around A 2 = 0.45 ± 0.1 GeV 2 . For the vector 
and the axial-vector quark-gluon-antiquark condensate we use a parametrization suggested 
in [15, 31]: 

(q{Qh»(-gA u (y))q(x)} = (y fl x u - g^iy ■ x))M 1 (x 2 ,y 2 ,(y - x) 2 ) 

+ {Vtf)u ~ 9^y 2 )M 2 (x 2 , y 2 , (y - x) 2 ) , 
{q{V>)lhlv.{-gA v {y))q{x)) = ze^yxMsix 2 , y 2 , (y - xf) 

with 

oo 

M,(x 2 , y 2 , z 2 ) -AiJJJda d/3 rf 7 fi(a, (3, 7) e H 2+/3 ^ 2 ) /4 . (A.2) 


Note that here the following abbreviation ^1,2,3 = A x (— 1,2,|) has been used. The 
minimal model of the nonlocal QCD vacuum relies upon the following ansatz: 

/, (a, 0, 7) = 6 (a - A) 6 ((3 - A) 6 (7 - A) (A.3) 

with A = A 2 /2. By construction, this type of model faces problems with the QCD equations 
of motion and the deficiency of the two-point correlator of the vector currents to satisfy the 
transversality condition. In order to fulfil the QCD equations of motion and at the same 
time to minimize the non-transversality of the VV correlator, an improved model of the 
QCD vacuum was suggested in [18] that uses a modified Gaussian ansatz, viz., 

jf P (a, (3, 7) = (1 + X t d x + Yidy + Yid z ) 6 (a - xA) 5 ((3 - yA) 5 (7 - zA) , (A.4a) 

where z = y, A = |A 2 and 

X l = +0.082 ; X 2 = -1.298 ; X 3 = +1.775 ; x = 0.788 , (A.4b) 
Fi = -2.243; Y 2 = -0.239 ; Y 3 = -3.166; y = 0.212. (A.4c) 

By virtue of the QCD equations of motion, these parameters satisfy the supplementary 
conditions: 

12 (X 2 + Y 2 ) -9 (X 1 + Y 1 ) = 1, x + y = l. (A.5) 

On the other hand, the vacuum condensates of the four-quark operators are reduced to a 
product of two scalar quark condensates (A.l) by virtue of the vacuum-dominance hypothesis 
[67]. 



21 



APPENDIX B: QCD SUM-RULE PARAMETERS 



The parameters to characterize the nonlocal condensates are A = A 2 /2 = 0.2 GeV 2 , 
(a s GG)/n = 0.012 GeV 4 and a s (qq) 2 = 1.83 • 10" 4 GeV 6 . The nonlocal gluon-condensate 
contribution §q(M 2 ) produces a very complicated expression. In analogy to the quark 
case, we model it by an exponential factor [13, 75]: $ G (M 2 ) = $q c (M 2 ) e " A 9 Q2/M4 with 
A 2 = 0.4 GeV 2 . 

In order to determine the best value of the threshold sq, we define a x 2 function for each 
value of Q 2 and sq : 



x 2 (Q 2 ,s ) 



N, 



M 



N M 

£ 

i=0 



£Q 2 F(Q 2 ,M 2 , So ) 



Q*F(Q 2 ,M 2 ,s ) 



\i=0 



2n 



;b.i) 



where we used M 2 = M 2 + i A Ml A M = (Ml - M 2 )/N M} N M = 20, and with e denoting 
the desired accuracy for x 2 — 1 (the actual value in the computation is e = 0.07 GeV 2 .) 



APPENDIX C: FACTORIZABLE PART OF THE PION FORM FACTOR IN 
FRACTIONAL ANALYTIC PERTURBATION THEORY 

Here, we provide some technical details about the FAPT analytization procedure in the 
computation of the factorizable part of the pion form factor using perturbation theory. In 
perturbative QCD at the NLO level and with the scale setting = Q 2 and /x| = const 
(where the subscripts F and R denote, respectively, the factorization and the renormalization 
scale), one obtains 



af\Q 2 ) 



n 2 



7T 



^ NLO (Q 2 ;/4) 



i^ pQCD ' (2) (Q 2 ;/4) = af\Q 2 )^°{q 2 ^l) 

where the superscript ^ in af (Q 2 ) denotes the two-loop order of the coupling, with 



(CI) 



Q 2 



LO/ ,,2\1 2 



[l + a^°(/i|) + a LO (/i| 



and 



^ NLO (Q 2 ;/4) = bort m (Q 2 ;ti) +^ FG \Q 2 ^l) + C F ^ F \Q 2 ^ 2 )ln 



L/4 



(C.2) 



. (C.3a) 



In Eq. (C.l), and below, and a\° are the LO Gegenbauer coefficients of the pion DA, 
whereas the individual contributions in Eq. (C.3a) read (with further details and explanations 
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given in [10, 52]) 



r 5 3 + (43/6)a£°(/4) + (136/15)a^°(^) 
3 + 



x 



l + 4°K)+a^(^) 
l + ^(/4) + a4°(/4)l S 



^• FG) (g 2 ;^) = {l5.67 + 4°(^) [21.52 - 6.22 4° (^)j 



.LO/..2 



(/4) [7.37 - 37.40 a£°(/4) - 33.61 a£°(/4 



Q 2 



2 r 
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^°(/4) 



i + 4 (/4) + ar(/4 



LO/,,2- 



(C.3b) 

(C.3c) 
(C.3d) 



The analyticized NNLO contribution (which includes radiative corrections of 0{a 2 s )) is then 
computed to be 



^ NLO (Q 2 ;^)-^ F )(Q 2 ;^) In 4 



■4 2) (g 2 ; 



A 2 



(C.4) 



Here A\{Q ), )> an d ^2 ; i(Q ) are the analytic images of ), a 2 (Q ), and 



« 2 (Q 2 )ln(g 2 /A^ CD ) [52, 70, 71]: 



A(Q 2 ) = A E [a s (Q 2 )] , 
A(g 2 ) e= A E [a 2 (Q 2 )] , 



£ 2; i(g 2 



«f(g 2 ) In 



g 2 



A 2 

QCD 



(C.5) 
(C.6) 

(C.7) 



with 



A E [f(Q 2 )} 



' Im/(-a) 
<y + Q 2 



da . 



{C.i 



This type of approach was analyzed in [71] and it was shown that the analytization procedure 
significantly reduces the dependence of the results on the factorization-scale setting. More 
recently [52], this procedure was further improved by setting /ip = Q 2 . This renders the 
logarithmic term to cancel out, but swamps the appearance of complicated expressions like 
A E [a p s (Q 2 )a2n (g 2 )a2°(g 2 )] with p = 1,2 and n,m = 0,1,2. Nevertheless, the obtained 
results are very close to those previously found in [71] and in [10] with = //p = Q 2 . The 
improved result F FAPT '^(g 2 ) was recently obtained in [52]. 
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